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Abstract 



Context. From a number of today known Post Common Envelopes Binaries (PCEB) only a handful has yet been observed 
at near-infrared (NIR) wavelengths and an even smaller number has modeled NIR light curves. At shorter wavelengths 
one has access to the cooler and larger components of these systems and has the chance to detect emission from its 
faint and heavily irradiated atmospheres. 

Aims. By modeling NIR light curves of PCEBs we intent to constrain their system parameters and study the properties 
of the system components. 

Methods. Here we present simultaneous NIR JHK S light curves of two PCEBs obtained with the 4m SOAR telescope. 
Results. KV Vel and T W Crv are long period (P or b = 8.6h and 7.9h, respectively) PCEBs with large irradiation effects. 
The results of light curve fitting provided solutions with inclination i = (47 ± 5)°, mass ratio q — 0.3 ± 0.1 and radius 
of the secondary R 2 /a = 0.24±o;°3 (where a is the orbital separation) for KV Vel, and i = (42 ± 9)°, q = 0.28 ± 0.04 
and 7?2/a = 0.22 ± 0.01 for TW Crv, respectively. For KV Vel, we obtain an average value for the albedo of the 
secondary star of a = 0.43, consistent in the J, H and iG-bands. For TW Crv, on the other hand, we obtain values of 
a j = (0.4 ± 0.1) and an = (0.3 ± 0.1) for the J- and //-bands, respectively. 

Key words. <keyword 1 - keyword 2 - keyword 3> 



1. Introduction 

Cataclysmic Variables (CVs) are binary systems where a 
low-mass dwarf (LMD) star overfill its Roche lobe and 
transfers matter to a more massive white dwarf (WD) at 
mass transfer rates typically of M ~ 10~ 10 — 10 -8 M Q i/r _1 . 
Since the orbital separation (and the Roche lobe) is ex- 
pected to increase when a lo wer mass sta r transfers matter 
to a more massive one (see iHellierl |2001| for an overview), 
the binary must continuously loose angular momentum to 
sustain mass transfer. CVs are thought to be the remnants 
of initially wide detached binaries where the more massive 
star evolves to a giant and expel its outer layers through 
a common-envelope (CE) configuration leaving a detached 
LMD+WD system. A detailed descrip t ion of this evolution- 
ary picture is given bv lTappert et al.l (|2007l ). 

Recently, a series of efforts has been made to unveil the 
evolutionary picture of CVs i.e, its angular-mo mentum loss 
histor y. Although there are alternative models, iKolb et all 
concluded that the disrupted magnetic braking 



model (DMB) is the most plausible explanation for CV 
evolution. Since the secondary star is responsible for the 
magnetic braking mechanism, understanding its properties 
is essential for the development of the DMB theory. By 
studying secondary stars of pre-CVs and PCEBs we are 
able to access stars in the same regime as those of CVs 
without the accretion complexity. 

In that sense we have collected near-infrared (NIR) 
JHK S light curves of a sample of PCEBs covering a range of 
orbital periods from around the period gap {P or b ~ 2 — 3h) 



up to P or i, ~ 8h, the largest orbital period suitable for full 
coverage within one night of (ground) observation. This pa- 
per reports the results of the analysis of the data from two 
of these systems. The observations and data reduction are 
presented in Sect. [5J Sample selection criteria and specific 
information about the objects of this study are given in 
Sect. El By modeling the light curve of these systems we 
are able to put constrains on their orbital system param- 
eters and, from the colors of each component, study their 
properties. In Sect. 0] we present the data analysis and dis- 
cussion of each target separately. Conclusions and further 
perspectives are presented in Sect. [S] 



2. Observations and reduction 

All observations were performed in queue mode with the 
4.1m SOAR Telescope at Cerro Pachon, Chi le, using the 
OSIR IS Infrared Imager and Spectrograph (jPogge et all 
1999) in imaging mode. A summary of the observations is 
presented in Table [TJ where t exp is the exposure time and 
ndith is the number of dithering positions. All runs were 
carried out under non-photometric conditions, with thin 
clouds and cirrus, resulting in variable sky transparency. 
Nevertheless, none of the observations had to be interrupted 
due to bad weather. 

Time-series of NIR photometry were performed quasi- 
simultaneously in J, H and K s , except in cases where 
the system was too faint and/or has too short orbital pe- 
riod that it was not possible to obtain good signal-to-noise 
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Table 1. Journal of observations. Date format is yyyy-mm-dd. 



Target 


Date of 


Bands 








ndith 


name 


observation 




J 


H 






KV Vel 


2007-03-24 


JHK a 


10s 


10s 


20s 


5 


TW Crv 


2007-04-24 


JH 


10s 


20s 




6 



data within the constraints given by the exposure time 
and the minimum time-resolution. The procedure of quasi- 
simultaneous photometry is basically to perform a set of 
acquisitions in one filter, change to the next filter, perform 
another another set and so on. This method warrants that 
any feature, lasting longer than an acquisition cycle, will 
be present in all light curves. 

In some cases, the required total exposure time to 
achieve the desired signal-to-noise ratio (S/N) saturates the 
profile of the target or of a field comparison star. In these 
cases we have to split a single exposure into more than one, 
and the final image is the sum of NC coadd images. For 
the targets discussed in this paper, the exposure time were 
short enough that no coadd procedure was necessary. 

Data reduction was performed using IRAlfl All frames 
were first trimmed and the non-linearity of the detector 
was co rrected by employin g the 3rd order polynomial solu- 
tion of iPogge et "all (|1999( 1 for the OSIRIS instrument. The 
frames were then divided by the normalized flat-field. A 
bad pixel mask was constructed from flat-field images and 
used to correct for bad pixels. Since the sky contribution to 
NIR light is considerably high, it is standard to perform a 
dithering procedure, nodding the telescope into ndith dif- 
ferent positions after each image is acquired. The sky level 
is obtained by taking the median of the ndith image. Each 
image is then iteratively shifted to match a reference image 
(usually the first image of the night) and aperture photom- 
etry is extracted for the target and all possible field stars 
(at least one star, aside the target, for each image field is re- 
quired) using a "variable aperture" (aperture radius = 2x 
FWHM). The filter was changed after each set of ndith 
acquisitions. 

Differential light curves (target star flux divided by com- 
parison star flux) were computed in order to account for sky 
transparency fluctuations during the night, and flux cali- 
brated from the 2MASS JHK S magnit udes of the field 
comp arison star and zero point constants (jSkrutskie et al.| . 
2006). 



3. Sample 

We selected a sample of southern PCEBs without accu- 
rate NIR light cur ve models. We searched in the catalog of 
iKube et alj f2002) for objects with confirmed orbital period 
of less than an observing night (P or b < 8h) and for systems 
where the exposure time required to achieve S/N = 50 
is less then Porb/50/N fi iters (with OSIRIS+SOAR, where 
N filters is the number of filters to be used for simultane- 
ous photometry) . The later is a requirement for light curve 
modeling. Since we are interested in secondary stars prop- 
erties "as-a-star" we do not exclude systems with subdwarf 



primaries. Some properties of the targets analyzed in this 
work are presented in Table [5] and discussed below. 

3.1. KV Vel 

KV Vel (or LSS 2018) is a non-eclipsing binary contain- 
ing a hot subdwarf primary (Ti = 77 000K) and a faint 
low-mass secondary star, and was the first central object 
of a planetary nebu la discovered to be a double-lined bi- 
nary (Drilling], |1985| ). iDrillind (|1985l ) provided ephemeris, 
radial velocity measurements and parameter determination 
of the system using UBV photoeletric measurements, IUE 
sp ectrophotometry and high resolution spectroscopy. 

IDrillind (|1985f) modeled the UBV light curves of KV 
Vel considering a simple heating model with spherical stars 
and no additional flux from the secondar y (bes ides that 
produced by heating). lLandolt fc Drillin g (1986) applied 
the same model to UBVRI light curves of KV Vel and 
pointed out that the quality of the fit degraded for wave- 
lengths longer than 5500A. Beyond this wavelength the 
secondary star contribution likely becomes significant and 
its distort ed shape could accou nt for the miss-fitting. In 
that sense. iHilditch et all (|1996l ) analysed the same data of 
lLandolt "fc Drilling] (Il986f ) with an improved model that ac- 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 



counted for the secondary flux and distorted shape as well 
as for illumination and atmospheric (e.g., limb- and gravity- 
darkening) effects. Although this improved model was able 
to provide good fit to the data, a problem wit h the modeling 
i llumi nation procedure was later reported by IHilditch et all 
(|2003f ). After inspection of the previous model, they found 
a missing factor of ir combined with a miss calculation of 
the surface normal vector of the emergent flux. These re- 
sulted in an underestimation, and possibly misshape, of the 
reflection effect, leading to unrealistic physical parameters 
from the fit. 

As stated by Landolt fc D rilling (1986J) and later by 
IHilditch et al.l (|1996h . modeling of infrared light curves of 
this binary would be key to our understanding of its prop- 
erties and of others of its kind, setting our motivation for 
this analysis. 

3.2. TW Crv 

TW Crv (or EC11575-1845) is a close binary similar to 
KV Vel. Both are long orbital period systems containing 
a hot subdwarf primary and a cool M dwarf secondary. 
A detailed descri ption of this system was presented by 
IChen et al.l {l995), who provided the first analysis of this 
close binary, combining UBVRIJHK photometry and opti- 
cal spectroscopy. Their light curves were used to derive an 
ephemeris for the system, 

T mid (E) = 2448661. 6049(±3) + 0.32762(±3) x E, (1) 

by fitting a sinusoidal curve to determine the time of max- 
imum light. A simple model of irra diated non-em itting 
spherical secondary, similar to that of Drilling (1985), was 
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Table 2. Properties of sample objects. Magnitudes were extracted from the 2MASS catalog. 



name 


Porb[d\ 


K 8 


J-H 


H — K a 


SP1 


SP2 


KV Vel 


0.36 (1) 




-0.116 l3) 


+0.644 (3) 


sdOB w 


M6±1V (4) 


TW Crv 


0.33 (2) 


13.3 (3) 


-0.022 (3) 


-0.016 (3) 


sdO< 2 > 


M?V (2) 



References. (1) Kilkenny et all (|1988h : (2) IChen et all ljl995l ): (3) ISkrutskie et all (|2006t ): (4) iDrillind lf l98 



used to model its optical light curves. As noted, this model 
produces poor fit to NIR light curves, since the assumption 
of a non-emitting secondary clearly fails at longer wave- 
lengths. 

Radial velocity curves of the system were constructed 
using the CIII/NIII blend at A4650A, providing velocity 
semi-amplitudes of Ki = (40 ±2) km/s and K 2 = (119 ±4) 
km/s for the prim ary and for the secon dary, respectively 
(jChen et al.Ul995D . As lChen et"al] (| 19951) pointed out, the 
measured K2 is actually the radial velocity of the irradiated 
face of the secondary, rather than the velocity of its cen- 
ter of light. They derived equations to correct that effect, 
as a function of the radius of the star, assuming uniform 
irr adiation. 

lExter et al.l ([2005) reported the analysis of high reso- 
lution (l.lA) time-resolved spectroscopy of TW Crv. They 
provided radial velocity measurements for both components 
of the system yielding values of Ki = (53 ± 2) km/s and 
K2 = (125 ±2) km/s for the p rimary and f o r the secondary, 
respectively. As in the case of IChen et al.l (jl995t ). the later 
is also the radial velocity of the centre of light and must cor- 
rected to the center of mass of the star before being used 
to derive binary parameters. 

For this analysis we have used new NIR light curves 
with higher S/N and lower exposure time and an improved 
model to analyze the data. 

4. Data analysis 

We analyzed the data wi th the aid of a light c urve m odel- 
ing code similar to the (|Wilson fc Devinnevl . 119711 here- 
after WD) algorithm. The code was developed by us 
to model the NIR light curves of detached and semi- 
de tached close binaries, acc ording to the routines outlined 
bv lKallrath fc Milond (1999). Our code was checke d against 
the W D code and also with the code developed bv IWatsonl 
IpOOl) . It was successfully applied to derive orbital param- 
eters of the CV IP Peg and, by extracting the light curve 
of the secondary, allowed the applic ation of eclipse map - 
ping techniques to its accretion disk (|Ribeiro et al.l . 12007). 
In the case of pre-CVs, the surface of the secondary is al- 
lowed to be any equipotential within its Roche lobe for a 
given filling factor, which gives the gravitational potential 
at the surface of the star with respect to the potential of the 
Roche lobe (// QC = ^2/^RLi > 1-0). Thus, a filling factor 
of unity corresponds to a Roche lobe filling secondary star, 
while ff ac > 1.0 represents under-filling components. 

The radiation field of the modeled component is modi- 
fied to account for gravity- and limb-darkening as well as for 
illumination effects. The latter considers a point-like source 
at the position of the primary star with bolometric irradi- 
at ion Ij r r . We have used the gravity-darkening coefficient 
of ISarnal (| 19891 ) for Roche lobe filling stars with convective 
envelopes ((3 = 0.05). Limb-darkening coefficients will be 
treated when discussing each target separately. 
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Figure 1. Orbital light curves of KV Vel, from top to bot- 
tom in J, H and K s bands, respectively. Model light curves 
are overploted as solid lines. The light curves are repeated 
in phase for better visualization. 



The parameters of the modeling procedure are: mass 
ratio (q = M2/M1), inclination (i), fluxes of the secondary 
(F2) and of the primary (Fi), irradiated intensity (Ii rr ) and 
filling factor (// ac )- Since we have no a priori information 
on the orbital separation and/or primary temperature with 
light curve modeling, we have used 7j rr to account for the 
amplitude of the irradiation effect (assuming an albedo of 
a = 1.0) and adopted the orbital separation as our distance 
scale. Once we know the amount of irradiation we can dis- 
cuss its origins. All light curves are modeled simultaneously. 
For a set of n light curves, q, i and ff ac are common pa- 
rameters while F2, F\ and Ii rr are particular to each data 
set. We proceed by minimizing the x 2 °f the model with 
respect to the data. The procedure works by first employ- 
ing a Simulated Annealing scheme to search for the region 
of best solut i on, an d then an amoeba minimization routine 
(jPress et al.l . Q.986) is employed to fine-tune the solution. 
Below we present and discuss the results of the modeling 
of the data for each target. 

4.1. KV Vel 

Our NIR light curves of KV Vel (Fig. [T|) are largely dom- 
inated by reflection ef f ects, and resembles its optical light 
curves (|Hilditch et all . Il996t ). The amplitude of the reflec- 
tion effect at NIR is AH ~ 0.7 mag against 0.55 mag in 
the V band. 

We fitted separate sinusoids to the JHK S data in order to 
measure the time of maximum light. The symmetric shape 
of the orbital hump with respect to the phase of maximum 
indicates that it is produced by uniform irradiation of the 
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Table 3. Times of maximum light and 0-C residuals for 
KV Vel. Our cycle timings is the mean of individual mea- 
surements of the JHK S bands. 



cycle 



T(max) 
(2400000H 



(0-C) a (O-C) 6 



Ref. 



-106 


45796.671 


+0.0084 


-0.0019 


1 





45834.52803 


+0.0106 


+0.0011 


2 


3 


45835.59908 


+0.0103 


+0.0009 


2 


829 


46130.5735 


+0.0032 


-0.0000 


2 


2834 


46846.5850 


-0.0118 


-0.0000 


3 


2845 


46850.5130 


-0.0122 


-0.0003 


3 


3018 


46912.2940 


-0.0130 


+0.0002 


3 


23382 


54184.53287 


+0.0045 


-0.0092 
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(a) With 
w W ith 
(1986). 



respect to the linear ephemeris of Eq. [5] 

respect to the linear ephemeris of Landolt & Drilling 



Ref erences. (1) [D rilling 
(3) iKilkennv et al.l i|1988 



CHI); (2) lLandolt fc Drilling] fl98rj ); 
(4) This work. 



secondary star by a source centered at the position of the 
primary star. In this case, we expect that the maximum of 
the light curve coincides with binary phase <p — 0.5, where 
the observer looks directly at the irradiated face of the sec- 
ondary star. Ho wever, if we phas e -fold the data according to 
the ephemeris of lKilkennv et all (|l988l ). the maximum light 
is displaced from its expected phase by Acf> ~ 0.01. The dis- 
crepancy can not be accounted for b y the uncertainty in the 
ephemeris of lKilkennv et all (|l988l ). 

Therefore, we combined our measured time of maximum 



with those of 


Drillina 


Kilkenny et al. 


(19881 



see Table [3]) to compute a revised 
ephemeris for KV Vel. The best fit least-squares revised 
linear ephemeris is, 

T max = HJD 2 445 834.5174(±4) + 0.3571205(±5) • E, (2) 

where E is the cycle number. The standard deviation of the 
data with respect to this ephemeris is a — 9.85 x 10 -3 , for 
a reduced chi-square of xt — 1-32. The (O-C) values with 
respect to this revised ephemeris are listed in Table [3] 

Although the revised linear ephemeris pr ovides a nice fit 
to our data, it yields a poor fit to the data of lKilkennv et all 
). These deviations prompted us to check whether a 



quadratic ephemeris provides a better fit to the whole d ata 
set. We applied the F-test proposed by iPringld (|1975| ) in 
order to determine the statistical significance of adding an 
additional term to the linear ephemeris. For this case we 
obtain F(l,6) = 2.8, with a statistical significance lower 
than 85% for the quadratic ephemeris. We conclude that 
presently there is no evidence of period changes in KV Vel. 

We have use d the previous determ ination of secondary 
star properties (jHilditch et all Il996|) to select square root 
limb-darkening coefficients from IClaretl (|l992h . The final 
set of parameters that best describe our data are shown 
in Table [?J and the corresponding light curves are shown 
as solid lines in Fig. [TJ Error estimation was performed 
through a Monte-Carlo simulation. 

Combining the inferr ed inclination with t he radial ve- 
locity measurements of iHilditch et al.l (|1996l ). we obtain 
masses of Mi = (O.7±O.2)M and M 2 = (O.3±O.1)M and 
an orbital separation of a = (2.1 ±0. 3)f? ( 7 1 . Our infe r red in - 
clination is slightly lower than that of Hilditc h et al.l |l996) , 
leading to higher masses for both components. The in- 



Table 4. Inferred parameters of KV Vel. 

Fluxes of the different components. 
J H 



K s 



F 2 (mJy) 0.40 ± 0.02 0.48 + 0.1 0.41 + 0.02 
Fi(mJy) 11.4 + 0.2 7.6 + 0.1 4.2 + 0.1 
hrr 26.6 + 0.2 15.6 + 0.1 10.1 + 0.1 
Simultaneously fitted parameters: 
i = (47 + 5)° q = (0.3 ± 0.1) f fac = (1.09 ± 0.02) a 

(o) The filling factor (f fac = $ 2 /$mi > 1.0) is the relative 
gravitational potential at the stellar surface. 



0.24ir o 5 3 or 



ferred filling factor yields a radiut0 of r2 
i?2 = 0-50tn n I fifT) fo r the secondary star. As noticed by 
Hildi tch et alJ (|l996l ). and underscored by our measure- 
ments, the secondary star of KV Vel is rather oversized 
for its mass in comparison to isolated main-sequence stars 
of the same mass — likely an indication that it is out of 
thermal equilibrium. 

The resulting mass for the primary star of KV Vel is only 
marginally consistent with the canonica l mass for subdwarf 
stars (M sdo = 0.4 7M^:lHan et al.ll2003D . The standard sce- 
nario proposed by lHan et al.ll2003 also predicts the forma- 
tion of stars with masses ranging from 0.3M Q to 0.8M Q . In 
this regard, the primary star of KV Vel may corre- 
spond to the rare case of the high mass end of the 
predicted subdwarf distribution. Nevertheless, we 
point out tha t the r adial vel ocity m easurements of 
KV Vel, from Hildi tch et al.l £[996), are potentially 
problematic. For instance, the lack of detailed mod- 
eling of the irradiation effect may conceivable hide 
important systematic errors. Therefore, further ra- 
dial velocity measurements, with proper modeling 
of irradiation, are required to solve this issue. 

We used the NIR colors of the low mass component of 
the binary to investigate its atmospheric properties and to 
constraint the distance to the system. The NIR colors for 
the secondary star of KV Vel are consistent with that of an 
M6V-M5V dwarf star and with those of a black body ra- 
diator with temperature Tbb = (3000 ± 100)K. Combining 
the colors and the inferred radius of the star, we infer a 
distance of d = (680 ± 60)pc to the system. In addition, 
we have also appli ed the stellar evolutionary models of 
iBaraffe et al.l (|1998| ). and consistently obtained a distance 
of d = (700 ± 100)pc to the system. The stellar models that 
provided the best-fit to the NIR flux distribution of the sec- 
ondary star are those of a star with a mass of M2 = O.2M 
or M 2 = 0.055M Q and ages of t = 10 6 5 yr or t = 10 8 5 yr, 
respectively. 

The resulting model properties of the stellar component 
are quite controversial. While and ag e of 10 85 yr is in closer 
agreem ent to that estimated by ISchreiber fc Gansickd 
(2003), it is quite difficult to reconcile the corresponding 
small mass of 0.055M Q with the dynamical solution of KV 
Vel. For example, the discrepancy between the mass and ra- 
dius of the secondary star would be even more pronounced 
and harder to explain. On the other hand, although a stel- 
lar component with 0.2M Q is in closer agreement with our 
dynamical solution, it is hard to reconcile the correspond- 
ing 10 6 5 yr age with the evolutionary stage of a PCEB. 
However, one might argue that, since the binary has re- 



2 here we list r2 = Rija, the radius of the secondary in units 
of the orbital separation 
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Table 5. Inferred parameters of TW Crv. 



Fluxes of the different components: 

J H 

F 2 (mJy) 2.4 ±0.6 2.1 ±0.4 

F d \mJy) 17.8 ±0.6 7.9 ± 0.3 

hrr 11.6 ±0.2 6.3 ± 0.4 

Simultaneouly fitted parameters: 
i = (41± 9)° g = (0.28 ± 0.04) f fac = (1.13 ± 0.01) 



cently evolved from a common envelope phase, the sec- 
ondary star is still out of thermal equilibrium (which is also 
indicated by the large radius discrepancy) and, therefore, 
its properties resemble those of young stellar atmospheres. 
We conclude that the M 2 = 0.2M© is the most plausible 
solution. 

4.2. TW Crv 

The light curves of TW Crv are also dominated by reflection 
effect, similar to those of KV Vel. Therefore, the analysis 
of its NIR light curve is similar to the analysis of KV Vel. 

We initially phase-folded the data using the ephemeris 
of Eq. [T] As in the case of KV Vel, the orbital minimum 
does not occur at phase (j) = as expected, but rather 
at (f> = —0.13 (earlier than expected). This is not surpris- 
ing given t h e rela tively low precision of the ephemeris of 
IChen et al.l (|1995l ) and the long time span between their 
and our observations. In analogy to the previous section, 
we fitted a sinusoid to the light curve to measure the time 
of minimum li g ht, an d combined the new timing with those 
of IChen et al.l (|l995h to revise the ephemeris of TW Crv. 

The best fit least-squares revised linear ephemeris for 
TW Crv is, 

T min (E) = 2448661. 6049 (±3) + 0.3276074(±2) • E. (3) 

The results of the light curve fitting analysis are listed 
in Table [5] and shown in Fig. [3J The filling factor leads to 
a radius of r 2 = 0.2 2 ± 0.01 for t he sec ondary star which, 
using equation (5) of IChen et al.l (|1995l ). gives a mass ratio 
of q = 0.30 ± 0.01, in good agreement with our light curve 
fitting value. 

Using eq. (4) of IChen et al.l |l995) together with out- 
values of i and r 2 we obtain Mt = (0.6 ± O.1)M for 
the total mass of the system, M\ = (0.5 ± 0.1)M Q and 
M 2 = (0.2±0.1)M Q for the stellar component masses, and 
an orbital separation of a = (1.7 ± 0.1)i?©. The same cal- 
culations performe d with the radial velocities provided by 
lExter et all (|2005h . result in a total mass of M T = (0.9 ± 
O.1)M , Mi = (O.66±O.O5)M and M 2 = (0.24 ± 0.05)M Q 
for the stellar component masses. In addition, we obtain an 
orbital separation of a = (1.9 ± 0.1)i2@ and a radius for 
the secondary star of R 2 = (O.4±O.1)R0. These results are 
obtained b y applying the sam e correction to K2 as in the 
case of the IChen et al.l (|l995jl results. 

Given the discrepa ncies between the resul ts using the 
radial velocity values of lChen et alJ (|l995h and lExter et al. 
(2005), it seems that the corrections de rived bvlCheiT et al] 
(1995) does not apply to the data of lExter et al.1 (|2005h . 
For instance, a non-negligible contribution from the un- 
de rlying secondary st ar to the narrow emission lines used 
bv lExter et al.l (|2005h invalidates the procedure. A detailed 
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Figure 2. Light curves of TW Crv in J and H with the cor- 
responding model fits (solid lines). The data were doubled 
in phase for better visualization. 



modeling of the impact of irradiation to the line shape may 
be necessary to clarify this issue. 

The radius of the secondary star obtained from the pre- 
vious analysis indicates that the secondary star is signif- 
icantly oversized for its mass. In addition, the flux ratio 
(F 2 (J)/F 2 (H)) of the secondary star is consistent with a 
black body of T — 4500K at a distance of d ~ 420 pc. 

5. Discussion and conclusions 

New NIR photometry of two long period PCEBs, KV Vel 
and TW Crv were presented and discussed. By measuring 
times of maximum light of KV Vel and TW Crv we were 
able to improve the system ephemerides. 

Figs.[T]and[2]shows our flux calibrated phase folded data 
together with corresponding model light curves. The model 
accounts for both the constant contribution of the primarjQ 
and of the illuminated distorted secondary, and also for 
inhomogeneities in the irradiation field of the secondary 
due to atmospheric effects. We were able to model the NIR 
light curves leading to well constrained system parameters 
(Tables H and [S]). 

From our light curve analysis it is also possible to in- 
vestigate the irradiation effect on the secondary surface. 
In contrast to the standard procedure usually adopted 
for albedo modelin g on these kind of binaries (p. ex., 
iHilditch et al.llT996l ). in which "known" physical parame- 
ters for the components are adopted and used to calculate 
the resulting light curve, our modeling procedure makes no 
a-priori assumptions on components properties. Yet, the 
use of previous measurements as initial set of parameters 
for light curve modeling is useful for accelerating the fitting 
procedure. Moreover, once we know the amount of irradi- 
ation with respect to the average surface intensity of the 
secondary (I s ) and, assuming an irradiating energy, we can 
estimate the albedo of the star. 

3 We note that the model cannot separate the constant con- 
tribution of the primary from that of any other source, like a 
shell or nebula, that may be contaminating the light from the 
system. 
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Table 6. Values for the albedo resulting from the light 
curve analysis. 



KV Vel 
TW Crv 



J 

0.43 ± 0.04 
0.6 ±0.1 



H 

0.44 ± 0.04 
0.5 ±0.1 



0.41 ±0.04 



For this analysis, the relation of I s with the albedo can 
be derived from, 



Rt 



1 + A t ^ 



F ' 



(4) 



where T{ is the modified effective temperature, TJ is the 
local temperature at the stellar surface, R t is the lo- 
cal reflection factor, A t is the bolometric albedo, F s is 
the incident flux and Ft is the undisturbed local flux 
(|Kallrath &; Milonel I1999D . Rearranging the terms in Eq. 
H]we obtain, 



T, 



Is -I 

F r 



where 



F s 

F — — 

r ~ Ft 



(5) 



We adopted ou r inferred radius of t he sec ondary star 
and t he values of iHilditch et al.l (|1996t) and IChen et al.l 
<|1995T> for the primary stars of KV Vel and TW Crv, respec- 
tively. We proceed by considering that the p rimary and sec - 
ondary stars in KV Vel have Ti = 77000K prillind . [l985h 
and T 2 = 3 400K (IHilditch et all 1 1996h . respectively, and 
that they irradiate as black bodied- With these assump- 
tions we are able to calculate the emergent (Ft) and inci- 
dent (F s ) flux at the surface of the star, and complete the 
calculations of Eq. [SJ The results are listed in Table [5] 

The resulting value for the albedo is strongly dependent 
on the assumed temperature for both stars. Therefore, in 
the case of TW Crv where the temperature of the primary 
is not well constrained, it is difficult to estimate the albedo 
reliably. The results in Table [5] are obtained adopting a 
temperature for th e primary of T\ = (105 ± 20) x 10 3 K 
(jExter et all l2005h and T 2 = 4 500K (Sect. g^J for the 
secondary. 

Although the final value for the albedo still depends 
on the adopted temperature and irradiation field of both 
components (i.e. whether or not it is valid to approximate 
their radiation by a black body) we were able to provide 
a good fit to the data by adopting standard values for 
gravity- and limb-da rkening, differently from the results of 
Hildi tch et alj (|1996[ ). Further investigation on system pa- 
rameters and details of the brightness distribution of the 
secondary are in demand for a better understanding of 
the system's properties and the physics of irradiated at- 
mospheres. 
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4 These assumptions are necessary in order to provide a way 
to measure the incident and emergent flux at the surface of the 
secondary. 
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